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Abstract 

Given a particular solution of a one-dimensional stationary Schrodinger 
equation this equation of second order can be reduced to a first order 
linear ordinary differential equation. This is done with the aid of an 
auxiliary Riccati differential equation. In the present work we show 
that the same fact is true in a multidimensional situation also. For 
simplicity we consider the case of two or three independent variables. 
One particular solution of the Schrodinger equation allows us to reduce 
this second order equation to a linear first order quaternionic differen- 
tial equation. As in one-dimensional case this is done with the aid of 
an auxiliary quaternionic Riccati equation. The resulting first order 
quaternionic equation is equivalent to the static Maxwell system and 
is closely related to the Dirac equation. In the case of two independent 
variables it is the well known Vekua equation from theory of pseudoan- 
alytic (or generalized analytic) functions. Nevertheless we show that 
even in this case it is very useful to consider not complex valued func- 
tions only, solutions of the Vekua equation but complete quaternionic 
functions. In this way the first order quaternionic equation represents 
two separate Vekua equations, one of which gives us solutions of the 
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Schrodinger equation and the other one can be considered as an auxil- 
iary equation of a simpler structure. Moreover for the auxiliary equa- 
tion we always have the corresponding Bers generating pair (F, G) , the 
base of the Bers theory of pseudoanalytic functions, and what is very 
important, the Bers derivatives of solutions of the auxiliary equation 
give us solutions of the main Vekua equation and as a consequence of 
the Schrodinger equation. Based on this fact we obtain an analogue of 
the Cauchy integral theorem for solutions of the Schrodinger equation. 
Other results from theory of pseudoanalytic functions can be written 
down for solutions of the Schrodinger equation. Moreover, for an am- 
ple class of potentials in the Schrodinger equation (which includes for 
instance all radial potentials), this new approach gives us a simple 
procedure allowing to obtain an infinite sequence of solutions of the 
Schrodinger equation from one known particular solution. 



1 Introduction 

Consider the one-dimensional static Schrodinger equation 

u" + vu = (1) 
and the associated Riccati equation 

y' + y^ = -V. (2) 

Equation (pQ) is related to the (O by the easily inverted substitution 

u' 

y = —- 

u 

Thus solutions of the Riccati equation Q are simply logarithmic derivatives 
of solutions of the Schrodinger equation and vice versa solutions of are 
logarithmic antiderivatives of solutions of ^ . The generalization of this fact 
for a multidimensional situation was obtained in |S] (see also |1U|). Among 
the peculiar properties of the Riccati equation stands out an important 
theorem of Euler, dating from 1760. If a particular solution i/q of the 
Riccati equation is known, the substitution y = yo + z reduces (j2j to a 
Bernoulli equation which in turn is reduced by the substitution z = ^ to a 
first order linear equation. Thus given a particular solution of the Riccati 
equation, it can be linearized and the general solution can be found in two 
integrations. As a consequence of this, given a particular solution of the 
Schrodinger equation the general solution can be found from a first order 
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linear equation. This can be seen immediately from the factorization of the 
one-dimensional Schrodinger operator 

d^ + v{x) = {d + yo{x)){d-yo{x)) (3) 

which is valid if and only if i/q is a solution of Q • 

In the present work we show that given a particular solution of a multi- 
dimensional stationary Schrodinger equation this equation of second order 
can be reduced to a first order linear quaternionic differential equation. For 
doing this we use a quaternionic factorization of the Schrodinger operator 
proposed in (see also ^7,) and the results on the quaternionic Riccati 

equation from (Hj and [IDj where it was shown that having a particular solu- 
tion of the quaternionic Riccati equation one can reduce it to a second order 
linear equation. Here we show that the similarity with the one-dimensional 
situation is much closer, and one particular solution is sufficient to reduce the 
quaternionic Riccati equation to a first order linear equation. The resulting 
first order quaternionic equation is equivalent to the static Maxwell system 
and is closely related to the Dirac equation. In the case of two independent 
variables it is the well known Vekua equation from theory of pseudoana- 
lytic (or generalized analytic) functions (see, e.g., [21, 0) EH, El)- We 
show that even in this case it is very useful to consider not only complex 
valued functions, solutions of the Vekua equation but complete quaternionic 
functions. In this way the first order quaternionic equation represents two 
separate Vekua equations, one of which gives us solutions of the Schrodinger 
equation and the other one can be considered as an auxiliary equation of 
a simpler structure. Moreover for the auxiliary equation we always have 
in explicit form the corresponding Bers generating pair (F, G) , the base of 
Bers' theory of pseudoanalytic functions, and what is very important, the 
Bers derivatives of solutions of the auxiliary equation give us solutions of 
the main Vekua equation and as a consequence of the Schrodinger equa- 
tion. Based on this fact, for example, we obtain an analogue of the Cauchy 
integral theorem for solutions of the Schrodinger equation. Other results 
from theory of pseudoanalytic functions can be written down for solutions 
of the Schrodinger equation. Moreover, for an ample class of potentials in 
the Schrodinger equation (which includes for instance all radial potentials), 
this new approach gives us a simple procedure allowing to obtain an infinite 
sequence of solutions of the Schrodinger equation from one known particular 
solution. 

Besides this introduction the paper contains four sections. In Section 2 
we introduce necessary notations from quaternionic analysis. In Section 3 we 
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prove a spatial generalization of the Euler theorem for the Riccati equation 
and show how a particular solution of the Schrodinger equation allows us to 
reduce it to a first order quaternionic equation. We observe that in the case 
of two independent variables this first order equation represents two separate 
Vekua equations. In order to apply theory of pseudoanalytic functions to 
the resulting Vekua equations, in Section 4 we introduce some necessary 
definitions and results from Bers' theory. Finally in Section 5 we show 
how all the machinery of this quite forgotten mathematical theory allows 
us to obtain surprising results for the Schrodinger equation starting with an 
analogue of the Cauchy integral theorem and including infinite sequences of 
solutions generated by one particular solution. 

2 Notations from quaternionic analysis 

We will consider the algebra EI(C) of complex quaternions or biquaternions 
which have the form q = qQ+ qii + q23 + gsk, where {qk} C C, and i, j, k 
are the quaternionic imaginary units. 

The vectorial representation of a complex quaternion will be used. Namely, 
each complex quaternion q is a sum of a scalar qo and of a vector q: 

q = Sc{q) + Vec(g) = go + q, 

where q = gii+g2j + 'J'3k. The purely vectorial complex quaternions (Sc((7) = 
0) are identified with vectors from C'^. Note that = — < q, q > where 
< •, • > denotes the usual scalar product. 

By we denote the operator of multiplication by a complex quaternion 
p from the right-hand side: M^q = q ■ p. 

More information on the structure of the algebra of complex quaternions 
can be found for example in fTIT or jl2j . 

Let q he a complex quaternion valued differentiable function of x = 
(xi,X2,X3). Denote 

„ . d . d ,5 

Dq = i^q+J-K—q + k^g. 

oxi 0x2 0x3 

This expression can be rewritten in vector form as follows 

Dq = — div q + grad go + rot q. 

That is, Sc{Dq) = — divq and Yec{Dq) = grad go + rot q. Let us notice that 
= -A. 

If go is a scalar function then Dqo coincides with grad go. The expression 
go/go will be called the logarithmic derivative of go- 
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3 Reduction of the Schrodinger equation to a first 
order quaternionic equation 

Consider the equation 

(_A + n)/ = (4) 

pfi Pi2. a2 

where A = ^ + ^ + / and u are complex valued functions. We 

assume that / is twice differentiable. Together with (jlj we introduce the 
following quaternionic equation 

Dq + q2 = _u (5) 

where q is a purely vectorial differentiable biquaternion valued function. 

Theorem 1 ^ For an arbitrary scalar twice differentiable function f the 
following equality holds 

{D + M'')iD-M'')f = i-A + u)f (6) 

if and only if h is a solution of |3p. 

Thus, given a particular solution of © the Schrodinger operator in (@J 
can be factorized. 

Theorem 2 Solutions of ^ are related to solutions of |3p in the follow- 
ing way. For any nonvanishing solution f of ^ its logarithmic derivative 

is a solution of ^ and any solution q of |3|j is a logarithmic derivative of 
the form ^ of a solution of 

Proof. A direct substitution into the equation Q shows us that for a 
nonvanishing solution / of Q) its logarithmic derivative is a solution of 
©. Now let us suppose that q is a solution of ©■ From the vector part of 
© we have that q is a gradient of some scalar function ^: q = grad^. Then 
q can be represented in the form ((T)) where / = e^. Substituting ((Tj) in © 
we obtain that / is a solution of Q). ■ 

Remark 3 Theorems^ and\^ show us that equation ^ is a generalization 
of the Riccati equation. We will call it quaternionic Riccati equation. 



5 



Lemma 4 '91 For a nonvanishing scalar dijjerentiable function e there ex- 
ists a one-to-one correspondence between solutions of the static Maxwell sys- 
tem 

div(eE) = 0, (8) 



rotE = 0, (9) 

and solutions of the equation 

{D + M^)F = (10) 

where 

Vector E is a solution of ^ if and only if the vector F = y^E is a 
solution of lll(J\) . 

Proof. The system (jSJ, © can be rewritten in the form 

^„ grade ^ 
DB =< ,E > . 

e 

Let us make a simple observation: the scalar product of two vectors p and 
q can be written as follows 

<P,q>=-^(PM + MP)q. 

Then we have 

1 grade 1 ^ (11) 

Note that 

1 grad e grad ^/e 

2 e ~ 

Then equation (|llj) can be rewritten in the following form 

-^D{^E) + Eh = (12) 

where h = D\/e/-y/e. Introducing the notation F = -y/eE and multiplying 
((T^ by ^/e we obtain the equivalence of the system 0-© to (fTn|l. ■ 
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Lemma 5 J^l Let h be an arbitrary particular solution of |3p (then as was 
mentioned above it is a gradient of some scalar function S^). The general 
solution of |3p has the form 

q = h + g, (13) 
where g = (gvad^)/^ and ^ is a general solution of the equation 

A^' + 2(grad^,grad^') =0, (14) 

or equivalently of 

div(e^^grad^') = 0. (15) 
Proof. Substituting (jO)) in © gives 

Z)g-2(h,g) +g2 = 0. (16) 
Note that the vector part of (|16|) is rot g = 0, so that 

g = grad <I> 

for some function <I>. If ^' = e*, this is equivalent to 

g = (grad^)/^'. 

Equation (|l(ij) . written in terms of is 

112 1 
-:p(grad^')2 _ -A^- - - (grad e, grad + —(grad^')^ = 0, 

so that (fro)) is equivalent to 

+ 2 (grad ^, grad = 0. 

Noting that 

div(e^^grad^') = 2e^^ (grad ^, grad ^')+e^^A^ = e^^(A^'+2 (grad ^, grad ^')), 
this equation can be rewritten in the form 

div(e^^ grad^*) = 0. 

■ 

Now we are ready to prove a generalization of the Euler theorem for the 
quaternionic Riccati equation. 



Theorem 6 (Euler's theorem for the quaternionic Riccati equation) Let 
h = grad^ be a particular solution of The general solution of the 

quaternionic Riccati equation has the form q = h + g where g = and ^ 
is obtained from the equation 

grad^' = e~^F (17) 
where F is the general solution of il(J\) . 

Proof. According to Lemma it is sufficient to prove tliat ^' is a solution 
of (|14|) (or what is the same of (|15|)) if and only if the vector F = grad ^ 
is a solution of ()1U|). Let us notice that if ^' is a solution of H14() then the 
vector E = grad ^' is a solution of the system 0-0 where e = e'^^ and vice 
versa if E is a solution of 0-® then it is a gradient of some function ^ 
which is necessarily a solution of (|14j) . Now due to Lemma |31 grad ^ is a 
solution of (jSJ-® if and only if F = grad^' is a solution of the equation 
(Tinll where h = ^ = grad^. ■ 

Thus, given a particular solution of the quaternionic Riccati equation, 
the general solution reduces to the linear first order equation (fTU]) . exactly 
as in one-dimensional situation. 

Using Theorem [21 we immediately arrive at the following result for the 
Schrodinger equation. 

Theorem 7 Let /o be a nonvanishing particular solution of ^ and F be 
the general solution of Iil0\) where h = Dfo/fo. Then the general solution 
f of ^ is the logarithmic antiderivative of q; = q, where q = h + g, 
g = and ^ is obtained from the equation 

F 

grad4' = --. (18) 
Jo 

Proof. From Theorem |21 we have that / is the logarithmic antiderivative 
of q (equation Q), where q is the general solution of (jSJ. For a particular 
solution /o of Q the vector 

h = — — = grad In /o 
Jo 

is a particular solution of Q. Due to Theorem [HI the general solution of Q 
has the form q = h-|-g, where g = and ^ is obtained from the equation 

F 

grad^' = e"'''-^"F = — 
Jo 

where F is the general solution of UlUj) . ■ 
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Remark 8 The logarithmic antiderivative of q always exists and can be 
obtained easily. Being a solution of |3p q is necessarily a gradient of some 
function $ which can be constructed analytically. Then f has the form 
f = Ce*, where C is a complex constant. 

Remark 9 From Theorem it follows that for any vector F, solution of 
lll(J\} with h = Dfo/fo, the vector F//o must be a gradient of some scalar 
function ^, and this is really true. Let us show that indeed rot(F//o) = 0. 
Note that this condition is equivalent to the equality Vec(i^(F//o)) = 0. 
Consider 

^(f) = >F-^F = -(F^ + ^F) = 2<F,^>. 
Jo JO Jo Jo Jo Jo 

Thus Vec(D(F//o)) = and hence the vector F//o is a gradient. 

Remark 10 Let us summarize the results of this section in the following 
chain of actions which one should follow in order to obtain solutions of ^ 
from solutions of mU\) . 

Given a nonvanishing particular solution /o of the Schrddinger equation 
^ we construct the vector h = Dfo/ fo and consider equation U(J\) . Taking 
a solution F of m(J\) we find ^ from m8\) . Then we construct the vectors 
g = D^/^ and q = h + g. q is necessarily a gradient of some scalar 
function <I>. Finding $ we finally obtain a solution of ^ as f = Ce*, where 
C is a complex constant. 

Thus given a particular solution of the Schrddinger equation the 
general solution reduces to the first order equation It is interesting to 

note that due to Lemma^ equation il(J\) is equivalent to the static Maxwell 
system. It is closely related (see ^ and 'llJ ) also to the Dirac equation as 
well as to the Beltrami fields which are solutions of the equation rot i+af = 
(see, e.g., /i/ and llSj). 

The following statement gives us the way to transform solutions of the 
Schrodinger equation into solutions of (|1U() . 

Proposition 11 Let fi be another nonvanishing solution of the Schrddinger 
equation Then the ratio = fi/fo is a solution of the equation 

div(/o2grad^') = (19) 

and the vector F = /o-D(/i//o) is a solution of equation mU\) where h = 
Dfo/fo. 



9 



Proof. By Theorem |21 the vector q = Dfi j f\ is a solution of © • Then 
consider the vector g from Lemma 13 

g = q- h= ^- :^ = Dih/fo)/ih/fo). 
Ji Jo 

Thus we have that g = D^/^ where ^ = /1//0, and by Lemma |SJ ^ 
satisfies (fT^ . 

From (UHl) we obtain that F = /o-D* = /o-D(/i//o) is a solution of (fTU)) . 

■ 

Let us consider equation ()1U|) for a biquaternion valued function p whose 
scalar part is not necessarily zero 

{D + M^)p = (20) 

and let us use the following representation for biquaternions: 

P = Pi + P2h 

where Pi = po + psk and P2 = P2 — Pik. Then D = Di + D2], where 
Di = 53k, D2 = — 9ik and h = Hi+H2^, where Hi = /13k, H2 = /i2 — /iik. 
Using these notations, equation ()20() can be rewritten as the following system 

DiPi - D2P2 + HiPi - H2P2 = 0, (21) 



D2P1 + D1P2 + H2P1 + H1P2 = 0. (22) 

Now let us suppose that both p and h do not depend on X3. Then system 
(|21jl . (|22|) turns into the pair of decoupled equations: 

D2P1 + IhPi = (23) 

and 

D2P2 + H2P2 = (24) 

which are nothing but Vekua's equations describing pseudoanalytic or gen- 
eralized analytic functions (see, e.g., jSj and [T^l. 

Remark 12 Here we should mention that in general the components pQ,...,p3 
as well as hi , /12 can he complex valued functions, hence Pi , P2 and H2 can be 
bicomplex. Nevertheless this detail is insignificant for what follows, because 
all results from Bers' theory which will be used in the subsequent sections are 
valid for bicomplex solutions also. Of course, when u in ^ is a real valued 
function we can consider real valued solutions of ^ only. In that case \2'J\) 
and \24\j are usual Vekua's equations. 
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Remark 13 In what follows we consider u and f in ^ being independent 
of X3. Then given a particular solution /o of the general solution reduces 
to equation \24\) (which in this case is equivalent to hlU\) ). Thus we are 
primarily interested in solutions of \24\) , and \2^J\) can be considered as an 
auxiliary equation. Nevertheless as we will see in Section 5 equations \23\) 
and \2Ji\l are closely related to each other. With the aid of Bers' theory 
solutions of \24^ can be obtained from solutions of I12'-J\) and it is interesting 
to notice that by construction we always have at least two solutions of \2S\) 
in explicit form. It is easy to see that the functions 

F=^ and G = /ok 

Jo 

are solutions of i2(J\) where h = Dfo/fo and consequently they are solutions 

ofmi- 



4 Some definitions and results from Bers' theory 

Bers' theory of pseudoanalytic functions was essentially developed in [Hj 
(see also (HI). It is based on the so called generating pair, a pair of complex 
functions F and G satisfying the inequality 

Im(FG) > (25) 

in some domain of interest O which may coincide with the whole complex 
plane. F and G are assumed to possess partial derivatives with respect 
to the real variables x and y. In this case the operators dz = ^ + 
and dz = — i-^ can be applied (usually these operators are introduced 
with the factor 1/2, nevertheless here it is somewhat more convenient to 
consider them without it) and the following characteristic coefficients of the 
pair {F, G) can be defined 

FGz - FzG , FGz - FzG 



FG-FG ' ' ' FG-FG 



_ FGz — FzG _ FGz — FzG 

where the subindex z or z means the application of or dz respectively. 

Every complex function w defined in a subdomain of 0, admits the unique 
representation w = <J)F + ipG where the functions (f) and ip are real valued. 
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Sometimes it is convenient to associate with the function w the function 
CO = + iip. The correspondence between w and uj is one-to-one. 

Bers introduces the notion of the (F, G)-derivative of a function w which 
exists and has the form 

w = (t)zF + il)^G = Wz- A(^F,G)'w - Bi^F^G)W (26) 

if and only if 

(hF + VrG = 0. (27) 
This last equation can be rewritten in the following form 

which we call the Vekua equation. Solutions of this equation are called 

[F, G)-pseudoanalytic functions. If w is (F, G)-pseudoanalytic, the associ- 
ated function uj is called (F, G)-pseudoanalytic of second kind. 

Remark 14 The functions F and G are {F, G) -pseudoanalytic, and F = 
G = 0. 

Definition 15 Let iF,G) and (Fi,Gi) - be two generating pairs in Q,. 
{Fi,Gi) is called successor of {F,G) and {F,G) is called predecessor of 
(FuGi) if 

«(Fi,Gi) = a(F,G) O-nd ^'(Fi,Gi) = -B[F,G)- 

The importance of this definition becomes obvious from the following 
statement. 

Theorem 16 Let w he an {F,G) -pseudoanalytic function and let (Fi,Gi) 
he a successor of {F, G). Then w is an {Fi, Gi) -pseudoanalytic function. 

Definition 17 Let {F, G) he a generating pair. Its adjoint generating pair 
{F,G)* = {F*,G*) is defined hy the formulas 

2F 2G 

F = = — ^— , G = —= — ^— . 

FG-FG FG-FG 

Theorem 18 

(F,G)** = (F,G), 

a(F*,G*) = -O-iFfi)-, ^(F*,G*) = -^(F,G)> 

i>(F%G*) = -B(F,G)> -B(F*,G*) = ~b{F,G)- 



12 



Lemma 19 is a successor of {F,G) then {F,G)* is a successor 

of{Fi,Gi)*. 

The (F, G) integral of w; on a rectifiable curve F is, by definition, 
j wd^^F^Q-jZ = Ke J F*wdz — iKe J G*wdz. 
Another important integral is also needed 



j G*wdz + iRe J 



'"^^(F,G)'2 = Re / G*wdz + iKe / F*wdz 



(we follow the notations of L. Bers). 

A continuous function w defined in a domain is called {F, G)-integrable 
if for every closed curve T situated in a simply connected subdomain of 0, 



Theorem 20 An {F,G)-derivative w of an {F, G)-pseudoanalytic function 
w is {F, G) -integrable and * f^^ wd{^p^Q-)Z = uj{zi) — uj{zq). 

The integral * J^^ wd^^p Q-^z is called (F, G)-antiderivative of w. 

Theorem 21 Let {F,G) be a predecessor of {Fi,Gi). A continuous func- 
tion is {Fi,Gi)-pseudoanalytic if and only if it is {F,G) -integrable. 



5 Applications of Bers' theory to the Schrodinger 
equation 

Let us return to equations (|23|) and (|24() which in a two dimensional case are 
equivalent to the quaternionic equation H20() . In order to use Bers' notations 
from the preceding section we rewrite ()23(l and H24|l in the following form 



w- = bw (28) 

and 

= bv (29) 

where z = x -\- iy, x = X2, y = xi and instead of the imaginary unit k we 
write i. It is easy to see that b = —H2 = —(hfo/fo ^-nd w = Pi, v = P2. 
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As was mentioned above f Remark [T^ for equation we know always 
two solutions 

F = ^ and G = ifo, (30) 
Jo 

which obviously fulfill (|25|) . Thus {F,G) is a generating pair corresponding 
to equation (j^Hj)- We have 

«{F,G) — U, O(^^c') — — 

4 - n R - ^^■^Q 

-4(F,G) - U, ^(F,G) 

According to Definition El the characteristic coefficients for a successor of 
(F, G) have the form 



«(Fi,Gi) - 0; ^Fi,Gi) 







/o 



Then due to Theorem 1161 if ii; is a solution of ()28() then its {F, G)-derivative 
is a solution of the equation 

= -bW, (31) 

but solutions of the last equation multiplied by i become solutions of (0^ 
and vice versa. Thus we obtain the following statement. 

Theorem 22 Let w be a solution of i28\) . Then the function 

. dzfo_. 

V = iw = lywz H 7 — w) 

Jo 

is a solution of \29{) . 

It is easy to see that according to Definition El 

F* = G* = fo 

Jo 



and 



b(F*,G*) = -B(F,G) = -b- 



Thus the {F, G) integral of a function W is defined as follows 

Wd,p^G)Z = - Re / -^Wdz - iRe [ foWdz 
Jr fo Jr 
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f W f 

= Im / —dz — iRe foWdz. 
Jr Jo Jr 

From Theorems 1201 and 1211 we obtain the fohowing result. 

Theorem 23 Let v be a solution of h2y\) in a domain Q. Then for every 
closed curve T situated in a simply connected subdomain of il,, 

Re / ^dz + ilm / fovdz = 0. (32) 
Jr Jo Jr 

Proof. For any solution v of H29|) the function W = iv is a solution of 
(|31j) . As pTj) corresponds to a successor of {F,G), by Theorem l2fll W is 
{F, G)-integrable. That is 

Im / ^dz -iRe [ foWdz = 0. 
Jr Jo Jr 

Now substituting iv instead of W we obtain . ■ 

In order to analyze the meaning of this result for solutions of the Schrodinger 
equation let us rewrite some statements from Section 3 in our "two-dimensional" 
notations. 

Consider the equation 

{-d,d^ + u)f = (33) 

where u and / depend on x, y, and z = x + iy. For simplicity we consider u 
and / being real valued functions. The corresponding Riccati equation (jSJ 
takes the form 

d^Q+\Q\^ = u 

where q = Qj. Equation (|T?H) turns into p9|) . F from ([T?Hl and v from (|29|) 
are related by the equality F = vj. Then Theorem |7| in a two-dimensional 
situation can be rewritten as follows. 

Theorem 24 Let /o be a nonvanishing particular solution of and v be 
the general solution of h2y\) where b = —dzfo/fo- Then the general solution f 
of iS.^} is obtained from the equation dzf = Qf, where Q = dzfa/ fo+dz'^ 
and ^ is obtained from the equation = "u/Zo- 

As was explained in Remarks 181101 given a solution v of H29|) . the corre- 
sponding solution / of can be constructed analytically. The procedure 
requires on two steps to reconstruct the potential function from its gradient. 

Prom Proposition ^2 we obtain the following statement. 
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Proposition 25 Let fi be another nonvanishing solution of Then the 

function 

v = fodMfo) (34) 
is a solution of \29^) . where b = —dzfo/fo- 

Having this precise relation between solutions of (jHS]) and (PH]) we are 
able to prove the following result. 

Theorem 26 (Cauchy's integral theorem for the Schrodinger equation) Let 
/o and fi be two arbitrary nonvanishing solutions of i'j'^) in a domain Q. 
Then for every closed curve T situated in a simply connected subdomain of 

n, 

Re / d,i^)dz + ilm I fld,{fy)dz = Q. (35) 
Jr /o Jt Jo 

Proof. Substitution of into gives us the result. ■ 

In ExamDle l31l we will give a nontrivial example illustrating this theorem. 
From Theorem!^ and Theorem 1261 we obtain an analogue of the Morera 

theorem for the Schrodinger equation (|33|) . 

Theorem 27 Let /o be a nonvanishing particular solution of The 
function fi is a solution of also if icl5\) is valid for every closed curve 
r situated in a simply connected subdomain of il.. 

Consider equation H27|) with the functions H3U() . It takes the form 

(h + ifo^z = 0. (36) 

Proposition 28 Let the function w = (j)F + 'ipG be {F, G)-pseudoanalytic 
corresponding to the functions i!^(J\) . that is 

w = ^ + ifoi' (37) 
70 

is a solution of i28\) . Then 

< grad</>,grad'(/' >= 0. 

Proof. The function (|37() is a solution of ()28|) if and only if H36j) is valid. 
Let us rewrite (|36|) in the form 

'Px - foi^y = (38) 
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4>y + = 0. (39) 

From (PU]) we have /q = —4>y/ipx- Substituting this expression into we 
obtain 

or 



In general, solution of or equivalently of the system (|5H)) . seems 
to be a difficult task. Nevertheless for a quite general class of functions /o 
we can obtain solutions of H38|) . H39() explicitly. Let us make the following 
suppositions. 



Condition 29 Let /q he a function of some variable p ■ fo = fo{p) such 
that is a function of p. We denote it by s{p) = 

The simplest example of such p is of course any harmonic function. An- 
other important example is p = \/ x'^ + y^- 

Consider the system (jSHJ, (jSH) and look for (/) being a function of p : 
(j) = 4>{p) (as we show below such solution always exists). Then 

= -^0', V. = %<P'. (40) 
Jo Jo 

For the solubility of this system we obtain the following condition 



which can be written as an ordinary differential equation 

,fo^ 

' f ' 

JO 



b" + is- 2$)6' = 0. 



From here we have 

^'{p) = e-'^^^f'oip) 

where S{p) = J s{p)dp. 

With the aid of H40|) we can reconstruct tp. Nevertheless we are interested 
nor in (j) neither in tp but in (p^ and "ipz instead. Having them we construct 



17 



the function v = i{(pzF + i^zG) = i{4>z/ fa + ^V'z/o) which gives us a solution 
of PH) . We have 

= = e-^/oV. (41) 

and 

V;, = = -^e-'pz. (42) 

/o 

Then we obtain the following solution of (|29|) 

z;i = i((/)^F + i)zG) = 2ifoe-^pz. 

In a much the same way we can construct another solution of ()29|) looking 
for Tp = ip{p). Then 

(f^x = foPyi^', 4>y = -foPxi^'- 

and ip' = e"^ / /q. Calculating (pz and we obtain 

(pz = ie-^Pz (43) 

and 

^z = -72rPz- (44) 

JO 

Thus we arrive at the following solution of (|29|) 

V2 = H — 7 ^ ^-prPzh) = -2e — . 

Jo Jo Jo 

Denote 

F/ = y and = ^. (45) 

Then lm{F jGi) = e~'^^ |gradp|^ > and hence we have a generating pair 
for H29|) in explicit form. Note that (Fj, Gj) is not a successor of {F, G) but 
a successor multiplied by i: {Fj, Gj) = i{Fi, Gi). 

It is interesting to see what are the new solutions /i and /2 of the 
Schrodinger equation corresponding to Ff and Gj. As the procedure 
of converting solutions of into solutions of requires on two steps 
to find a potential function by its gradient, it will be more illustrative to 
consider here an example. 
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Example 30 Let u{x,y) = + y^. Then a particular solution of \8i^) can 
he chosen in the form 

fo{x,y)=e^y. 

Obviously p{x, y) = xy being a harmonic function satisfies Condition \2!A 
Then Fj = ie^y{y — ix) and Gj = —e~^y{y — ix). Returning to the notations 
of Section 3 (see the beginning of Section 5) we have Fj = he'^^^'^ {xi—k.X2) = 
e^^^^{x2 + kxi) and Gj = — e~^^^^(xi — kx2). Then the corresponding pair 
of solutions of mU\) has the form 

Fi = Fa = e^i^H-xii + X2j) and Fa = G/j = -e-^'i^^xai + xij). 

Next step (see Remark \1(J\) consists in finding the corresponding functions 
^1 and from equation I118\) . Thus we should reconstruct ^'i and ^'2 from 
the equalities 

grad ^I'l = — xii + X2j and grad ^'2 = — e~^^^^^(x2i + xij). 

Using the standard formula for finding the potential function from its gradi- 
ent we obtain 

^1 = -^(x? - - Gi) and ^2 = \{e-^^^^^ + G2) 

where Gi and G2 are arbitrary constants. 

Now we can construct the vectors gi = grad and g2 = grad ^'2/^'2-' 

2 G 

Noting that h =-D/o//o = X2i + a;ij we obtain two solutions qi = h + gi and 
q2 = h + g2 for 

2xi 2x2 



qi = (X2 + 2 7r)i + (^1 



x^ X2 Gi x^ X2 G\ 

and 

2C2X2 , , 2C2X1 

Now we find the functions and ^2 which are solutions of the equations 
grad<I>i,2 = qi,2- 

$1 = In \x1 — X2 — Ci| + X1X2 + C3 
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and 

$2 = xiX2 + In |e-2^i^2 + + C4. 
T/ien the corresponding solutions fi and /2 of have the form 
A = e*^ = d,{xl -xl- C7i)e^i^^ = _ - COe^^' 

and 

/2 = = d2(e-"^"^ + C2e"^"2) = ^2(6""^ + ^26^^') (46) 

where di and d2 are arbitrary constant. 

Thus starting with a particular solution of it ^5)) we constructed two classes 
of solutions for the same Schrddinger equation. 

Example 31 In order to illustrate the Cauchy integral theorem for the 
Schrddinger equation let us use the following two particular solutions from 
the preceding example. Let fQ{x,y) = e^^ and as f\ we choose the function 
from i46{ ) when C2 = and d2 = 1, fi{x,y) = e~^^. Both fo and fi are 
solutions of with the same potential u in a whole plane. Thus we can 
apply Theorem A2fA and consider T being for example a unitary circumference 
with centre at the origin. Then 

Re I d,{fy)dz + ilia [ fld,{fy)dz 
Jv JO Jt jo 

= 2ReJ (-y + ix)e"^^^dz + 2ilm j {-y + ix)dz 

/'27r /'27r 

= 2Re/ i(-sinr + zcosr)e"^™''''''™^dr + 2iIm / i(- sinr + i cos r)(ir 
Jo Jo 

-2 / cosT-e-^^°''^'''"^^dT-2i f smrdr. 



It is easy to see that both integrals are really equal to zero. 

Let us calculate the characteristic coefficients for the pair (Fj,Gi) de- 
fined by (Uni)- We have 



a(F,,Gi) = 0' hFi,Gi) 



dzf 







fi 







I „2 



. Pzz Pzz „ foPz 
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Consider the equation 

i-Fj + ip^Gi = 

1 1 

where (p and -0 are real valued functions. It has the form 

iff^-i;- = (47) 

or as a system 

k + foiy = (48) 

iy - fok = 0. (49) 

Comparing with the system (|38|) . H39() we note that solutions of that system 
can be transformed into solutions of the system (|48|). (|49|) in the following 
way 

1 1 

(j) = il) and V = ~<?^- 



Let us calculate and using Definition El 

F; = -:^°^ and '^"^ 



P2 /O/Oz 

Consider the equation 

2 2 

where (f) and are real valued functions. It has the form 

iff(t^-i^ = ^. (50) 

Observe that it coincides with (|47l) . Thus we obtain that the function v = 
(pFj + ipGi is (F/, G/)-pseudoanalytic iff the function W = (j)Fj + 'tpG*j is 
{Fj , G^)-pseudoanalytic. 
Let us calculate 

That is the characteristic coefficient 6 of a successor of {Fj,GJ) is equal to 
—9zfo/fo = ^{F,G)- Thus, {F,G) is a successor of {Fj,G*j). This important 
observation opens the way to obtain an infinite set of solutions of the original 
Schrodinger equation (|33|) if /o satisfies condition 1291 Namely, we start with 
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a solution of H28|) . for instance with F. Then its {Fj , G|)-antiderivative gives 
us an (Fj*, G})-pseudoanalytic function, more precisely the corresponding 
functions (f) and ip such that the function W = (pFj + ^l^GJ is {Fj,GJ)- 
pseudoanalytic. Then we take these real valued functions cp, ■0 and consider 
the function v = cpFj + which is (F/, G/)-pseudoanalytic. That is v 
satisfies ()29() and hence the vector F = vj with the aid of the chain of 
actions described in Eemark llOl can be transformed into a solution of . 
Taking the {F, (j')-antiderivative of iv we obtain another solution of ()28() and 
can start this cycle again. We can represent schematically this procedure 
for obtaining an infinite sequence of solutions of H29() and consequently of 
(|.S3j) as the following diagram. 

Solutions of dSHl) 

/ 

(F,G) ^ iFj,Gj) 

\ T 

iF*,G}) 

Let us consider how this procedure works on the following example. 

Example 32 Here we use the same u and /o as in Example Then 
F = I//0. We have that = -6^2^/(7/ - ix), G} = -ie-'^y/iy - ix). 
Consider 

* / Fdi^p*^G*~^z = ^Q \ FGjdz + iRe / FFjdz 
Jo ' Jo Jo 

= - Re [ e"^""^*' {yt - ixt) {x + iy)dt + z Re [ (xt + iyt) (x + iy)dt 
Jo Jo 



-2xy . (X^ - y' 



2^ 



2 2 

Thus we have (p = (e~^^^ — l)/2 and = (x^ — y'^)/2. It is easy to check 
that in fact (pzF^ + ipzGJ = 0, so the function 

W = 0F; + i;G*i = V , (e-"^ - e^^ + ie-^^^ix^ - y^)) 

2[y — IX) ' 

is {Fj, G*j)-pseudoanalytic. Now we can construct an (Fj, Gj)-pseudoanalytic 
function as follows 

v = (t>Fj+ ijGi = (e-^f (x^ - y2) _ i(e-^s/ _ g^^/)) . 
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Thus V is a solution of h2y\) and applying the procedure described above it 
can be transformed into a solution of h'J'Jd . 

Now multiplying the function v by i we obtain a solution of the 
{F, G)-antiderivative of which gives us a new solution of and the cycle 
starts again. 

Thus starting with a particular solution F of i28\) we obtain an infinite 
sequence of solutions of iV.V)) . 
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